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ABSTRACT 

We present new observations of a transit of the 1 1 1.4-day-period exoplanet HD 80606b. Due to this long orbital period and to the 
orientation of the eccentric orbit (e = 0.9), the HD 80606b's transits last for about 12 hours. This makes practically impossible the 
observation of a full transit from a given ground-based observatory. Using the Spitzer Space Telescope and its IRAC camera on 
the post-cryogenic mission, we performed a 19-hour-long photometric observation of HD 80606 that covers the full transit of 13-14 
January 2010 as well as off-transit references immediately before and after the event. We complement this photometric data by new 
spectroscopic observations that we simultaneously performed with SOPHIE at Haute-Provence Observatory. This provides radial 
velocity measurements of the first half of the transit that was previously uncovered with spectroscopy. This new data set allows the 
parameters of this singular planetary system to be significantly refined. We obtained a planet- to-star radius ratio R p /R, = 0.1001 ± 
0.0006 that is more accurate but slightly lower than the one measured from previous ground observations in the optical. We found 
no astrophysical interpretations able to explain such a difference between optical and infrared radii; we rather favor underestimated 
systematic uncertainties, maybe in the ground-based composite light curve. We detected a feature in the Spitzer light curve that could 
be due to a stellar spot. We also found a transit timing about 20 minutes earlier than the ephemeris prediction; this could be caused 
by actual transit timing variations due to an additional body in the system, or again by underestimated systematic uncertainties. The 
actual angle between the spin-axis of HD 80606 and the normal to the planetary orbital plane is found to be near 40° thanks to the 
fit of the Rossiter-McLaughlin anomaly, with a sky-projected value A = 42° ± 8°. This allows scenarios with aligned spin-orbit to 
be definitively rejected. Over the twenty planetary systems with measured spin-orbit angles, a few of them are misaligned; this is 
probably the signature of two different evolution scenarios for misaligned and aligned systems, depending if they experienced or not 
gravitational interaction with a third body. As in the case of HD 80606, most of the planetary systems including a massive planet are 
tilted; this could be the signature of a separate evolution scenario for massive planets in comparison with Jupiter-mass planets. 

Key words. Planetary systems - Techniques: radial velocities - Techniques: photometry - Stars: individual: HD 80606 



1. Introduction 

Among the more than 400 extrasolar planets that have been 
found so far, the giant planet orbiting the G5 star HD 80606 
is certainly a unique case. Its eccentricity is particularly high: 
e = 0.93. Only one known planet possibly has a higher eccen- 
tricity, namely HD 20782b; however its high e-value is still to 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, and with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. 



be confirmed as it stands on one measurement only (O' Toole 
et al. 120091 . The comet-like orbit of HD 80606 b was weU es- 
tablished in its discovery paper by Naef et al. (2001) and has 
been largely confirmed by subsequent observations. Together 
with its 111 .4-day period, the high eccentricity of HD 80606b 
put it on an extreme orbit: during its revolution, the planet expe- 
riences the strongly irradiated regime of a "hot Jupiter" at peri- 
astron (0.03 AU), and milder conditions at apoastron (0.87 AU), 
around which it spends most of its time, and where the planet ap- 
proaches the inner boundary of the habitable zone (see Fig . [TJ . 

As long as the inclination i of the orbit with the line of sight 
was unknown, those parameters originally implied a probabil- 
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Fig. 1. Schematic view of the HD 80606 system. The red dots 
show the positions of the planet each 24 hours. 



ity of 1% for the planet to be transiting. In spite of this tenuous 
probability, an amazing chance makes the planet HD 80606b ac- 
tually transits its parent star, as seen from Earth, each 111.4 days. 
This is particularly advantageous as numerous crucial studies 
can be performed using photometry or spectroscopy when a 
planet passes in front its parent star (planetary transits) or behind 
it (planetary eclipses), especially in this case where the host star 
is bright (V = 9.1) and nearby (d = 60 pc). The fortunate tran- 
siting nature of HD 80606b was established in February 2009 
from the detection of a transit reported from ground obser- 
vations, independently by Moutou et al. (120091) from pho- 
tometric and spectroscopic data, and by Garcia-Melendo & 
McCullough (2009) and Fossey et al. ( 120091 ) from photomet- 
ric measurements. Additional observations of transits were later 
reported by Winn et al. d2009aD and Hidas et al. d20TDb . The 
February 2009 observations followed the planetary eclipse dis- 
covery reported a few months earlier by Laughlin et al. (2009) 
from Spitzer photometric observations at 8 fim during a 30-hour 
interval around the periastron. Among the known transiting plan- 
ets, HD 80606b has the longest period and the most eccentric 
orbit. The second most extreme transiting planet is HD 17156b 
(P = 21.2 days and e = 0.67). Furthermore, HD 80606b is also 
the most distant planet from its parent stars when it transits: 
0.3 AU, against 0.05 AU or less for all other known casefl 

In addition to the photometric detection of the HD 80606b's 
transit, Moutou et al. ((2009) presented its spectroscopic detec- 
tion using the Rossiter-McLaughlin effect, from radial veloci- 
ties measured with the SOPHIE spectrograph at Haute-Provence 
Observatory (OHP), France. The Rossiter-McLaughlin effect is 
an apparent distortion of the stellar lines profile due to the tran- 
sit of the planet in front of the rotating star. From the SOPHIE 
measurements, Moutou et al. (2009) have shown the first ev- 
idence for a spin-orbit misalignment, i.e. the orbital plane of 
the planet HD 80606b is not perpendicular to the spin-axis of 
its host star. Using additional photometric data of the February 
2009 event allowing a better constraint on the transit duration 
together with a combined analysis of the whole data set, Pont et 
al. d2009) refined the parameters of the system. They confirmed 
the spin-orbit misalignment from the Rossiter-McLaughlin dis- 
tortion detected with SOPHIE and provided a measurement of 
the sky-projected angle between the planetary orbital axis and 
the stellar rotation axis: A ~ 50°, with the confidence interval 
[14° - 111 ] - see also Gillon (I2009al ). Thanks to new photo- 
metric and spectroscopic observations of the June 2009 transit, 



1 See also the 95-day-period transiting exoplanet CoRoT-9b an- 
nounced after the submission of the present paper (Deeg et al. l2010t . 



Winn et al. (2009a) thereafter reduced even more the confidence 
interval to [32° - 87°]. Thus, the spin-orbit misalignment of the 
HD 80606 system is now well established. HD 80606 is the com- 
ponent of a wide binary with HD 80607; the projected separation 
of the system is about 1000 AU. The peculiar orbit of HD 80606b 
could thus result from Kozai mechanism and tidal dissipation 
(see, e.g., Wu & Murray 2003), which can pump the eccentricity 
and the inclination. 

Due to the long orbital period and to the orientation of the 
eccentric orbit, the duration of the transit of HD 80606b is about 
12-hour long. This should be compared to the transit duration 
of other known transiting exoplanets, which typically lasts less 
than five hours. The transit duration of HD 80606b is even longer 
than transits of Mercury or Venus through the Solar disk as seen 
from the Earth. It is thus practically impossible that a full transit 
of HD 80606b matches the duration of an observation night from 
ground. In addition, data secured before and after the transit are 
mandatory to obtain an accurate transit light curve, so the full se- 
quence for HD 80606 lasts longer than a night for ground obser- 
vations. Observing an entire transit of this exoplanet is thereby 
challenging and only portions of a transit could be observed from 
a given ground-based telescope. This was the case of all the ob- 
servation campaigns reported above which covered only frac- 
tions of transits. These fragmented observations induce signifi- 
cant uncertainties in the parameters derived from their fit. 

We present here the first full photometric observation of a 
transit of HD 80606b. We secured it on 13-14 January 2010 with 
the Spitzer space observatory using the IRAC camera at 4.5 pm 
in post-cryogenic mission. Thanks to its Earth-trailing heliocen- 
tric orbit (Scoupe et al. 2006), Spitzer allowed us to continuously 
observe during 19 hours, enabling the coverage of the whole 
12-hour-long transit, as well as off-transit references immedi- 
ately before and after the event. We complement this photomet- 
ric data by new spectroscopic observations that we simultane- 
ously performed with SOPHIE at OHP. This provides radial ve- 
locity measurements of the first half of the transit, a part that was 
up to now uncovered with spectroscopy. Indeed, observing the 
full 12-hour-long transit is even more difficult in spectroscopy 
than in photometry, as the amplitude of the Rossiter-McLaughlin 
for HD 80606 is about 10ms 1 whereas Northern instruments 
able to measure radial velocities with the required accuracy are 
sparse. We also performed a ground-based photometric monitor- 
ing of HD 80606 during January 2010. All together, the data of 
this observational campaign allows the parameters of this plane- 
tary system to be additionally refined. 

The observations and data reduction are presented in Sect. [2] 
and [3] for Spitzer and SOPHIE respectively. The ground-based 
photometry in presented in Sect. |4] The analyses and the re- 
sults are presented in Sect. [5] before discussion and conclusion 
in Sect.|6]and [7] 

2. Spitzer photometry 

2.1. Observations 

We obtained Spitzer Director's Discretionary Time (DDT pro- 
gram #540) to observe the January 2010 transit of HD 80606b. 
This transit was the first observable with Spitzer since the dis- 
covery of the transiting nature of HD 80606b in February 2009. 
As Spitzer has exhausted its cryogen of liquid coolant on 15 May 
2009, our observations were performed during the first months 
of the Spitzer's warm mission. Only the two first infrared chan- 
nels of the IRAC camera (Fazio et al. 2004) are available in the 
post-cryogenic Spitzer. They are centered at 3.6 and 4.5 pm and 
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cannot be observed simultaneously. We chose to observe only in 
one of the two channels in order to avoid repointing the telescope 
during the transit. This reduces overheads and ensures that the 
target remains on the same part of the detector during the full ob- 
servation sequence. We opted for channel 2 at 4.5 pm since it has 
the lowest noise properties. This wavelength also has a smaller 
limb-darkening effect for the star. We did not use pointing dither- 
ing, here again to maintain the target as much as possible on the 
same location on the detector in order to reduce systematic ef- 
fects due to imperfect flat-field corrections and intra-pixel sensi- 
tivity variations. 

We used the Subarray mode of IRAC which is well adapted 
for bright targets. Only a 32 x 32-pixel part of the detector is used 
in this mode; this covers a small 38 x 38 arcsec 2 field of view 
(pixel size of 1.2 arcsec), compared to the IRAC Stellar mode 
that uses the full 256 x 256-pixel field. As the stellar compan- 
ion HD 80607 is located only 17 arcsec to the East of HD 80606, 
putting the two targets on this small field of view would have 
imply that their point spread functions (PSF) fall near the edges 
of the detector. Such configuration is risky for accurate photom- 
etry. We preferred to let HD 80607 off the field of view in order 
to maintain HD 80606 on the detector, far from its edges. So we 
chose to put our target at the default pointing position in the cen- 
ter of the Subarray field of view. This position is not on nor right 
next to any known hot pixels. 

The observations were secured between 2010 January 13 at 
18 h and January 14 at 13 h (UT). We acquired 2488 consecutive 
images during a total of 19 hours. Each image was split up into 
64 frames of 0.36 second exposure each, taken back-to-back. 
We obtained a total of 159232 frames during 15.9-hours effec- 
tive integration time with an overhead of 2-second between each 
images (15 % overheads in total). Such high efficiency is reached 
thanks to the use of the Subarray mode and could not be achieved 
in Stellar mode for bright targets requiring short exposure times. 
With a flux density of ~ 200 mJy at 4.5 pm for HD 80606, frame 
exposure time of 0.36 sec clearly avoids saturation of the pixels. 
The intensities recorded in the brightest central pixel are around 
10 000 electrons. 

2.2. Data reduction 

We used the BCD files (Basic Calibrated Data) of each of 
the 159232 frames as they are produced by the Spitzer/IRAC 
pipeline. It includes corrections for dark current, flat fielding, 
pixels non-linearity, and conversion to flux units. To extract the 
light curve, we used tools and methods we developed in Desert 
et al. (120091 [2Q10l >. We find the center of the PSF of the star to 
a precision of 0.01 pixel using the DAOPHOT-type Photometry 
Procedures, GCNTRD, from the IDL Astronomy Librar)Q, which 
computes the stellar centroid by Gaussian fitting. We used the 
APER routine to extract the raw flux of HD 80606 on each frame 
from the computation of a simple aperture photometry using a 
radius of 4.0 pixels, which optimize the signal-to-noise ratio of 
the transit light curve. The flux integrated on these 50 pixels is 
~ 52 500 electrons. It has been corrected from the background 
level of 14.40 ± 0.05 electrons/pixel that has been estimated 
from a sky annulus with radii of 9 to 12 pixels. The centroid 
of HD 80607 is located outside the field of view but a small con- 
tribution of its flux is detected on an edge of the detector. When 
estimating the background level, we took care to use only pix- 
els where the HD 80607's flux is negligible. It is also negligible 
by comparison to the background level in the 4-pixel radius we 



used for the HD 80606 photometry. Finally, the uncertainty in 
the background level is negligible for all the results presented in 
this paper. 

The Spitzer/IRAC photometry is known to be systematically 
affected by the so-called pixel-phase effect. This effect produces 
an oscillation of the measured raw light curve due to the Spitzer 
telescope jitter and the intra-pixel sensitivity variations on the 
IRAC detector (see, e.g., Charbonneau et al. [2005, Reach et 
al. 120061 Mo rales-Calder on et al. l2006l Ehrenreich et al. l2007l 
Desert et al. 2009, 2010). Measurements of the centroid posi- 
tion of the target on the detector and its variations could be used 
to de-correlate the pixel-phase effect on the light curve. We use 
here the method presented in Desert et al. (2009); it has the form 
^coir = F[l+K l (x-x ) + K2(x-x ) 2 + K 3 (y-y ) + K 4 (y-y ) 2 + 
K$(x - xo)(y - yo)], where F and F con - are the fluxes of the star 
before and after the pixel-phase effect correction, and (x - xq) 
and {y - yo) are the position in pixel of the source centroid on 
the detector with respect to the pixel pointing position, located 
on [^o^yo]- Our determination of the centroid position shows a 
+0. 05-pixel oscillation with a period of ~ 70 minutes and a lin- 
ear drift during the 19-hour sequence of 0.1 and 0.2 pixel in the 
x- and y-direction, respectively. 




10.4 10.6 10. i 

HJD - 2455200 



Fig. 2. Spitzer photometry of HD 80606. From top to bottom, the 
three panels show respectively the raw photometry, the photom- 
etry after correction of the pixel-phase and ramp effects, and the 
residuals. The black points are the measurements (binned per 5) 
and the red solid lines are the fits. The vertical dotted lines show 
the mid-transit, the vertical dot-dashed lines show the first and 
fourth contacts, and the vertical dot-dot-dashed lines show the 
second and third contacts. 
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In addition to the 64 frames of the image # 873 that are cor- 
rupted, we iteratively selected and trimmed 1037 outliers by 
comparison to the fit of the light curve with a transit model. 
Frames were considered as outliers when they were above 10 
to 3 <x, this value being reduced by 0. 1-cr steps at each iterations. 
We binned the obtained transit light curve by a factor of five in 
order to obtain a better computing efficiency without losing in- 
formation for the pixel phase effect. Most of the results presented 
below were obtained using the binned transit light curve. 

We tested several radii for the aperture photometry, sev- 
eral areas for the background measurement, several procedures 
for the centroid determination and the outliers rejections. The 
adopted procedure reported above is the one producing the 
smallest errors, but all of them produced similar results. 

The upper panel of Fig. [2] shows the raw Spitzer light curve 
of HD 80606 after this extraction. It clearly shows flux varia- 
tions with a period of ~ 70 minutes that are due to the pixel- 
phase effect and telescope jitter at this period. Its peak-to-peak 
amplitude represents ~ 1 % of the flux, which is larger than 
the effect seen in the channel 2 of IRAC during the cryogenic 
Spitzer. Additional variations at higher frequency and with a 
lower amplitude than the one of the 70-minute oscillation are 
seen as well, as we fit each individual frame. They could be due 
to short-term jitter of the satellite which apparently are not pe- 
riodic. Filtering these high-frequency variations does not signif- 
icantly change our results. Finally, a slope in the out-of-transit 
baseline is also seen. Such detector ramp was observed in chan- 
nels 3 and 4 in pre-cryogenic Spitzer but was not as significant 
in channel 2 at that time (see, e.g., Deming et al. 2005 Knutson 
et al. 120071 Desert et al. 120091 l2"0"T0l ). We discarded the first 
thousands of frames that were the most affected by the ramp 
effect (we tested several limits; see Sect. 15.1. 41 then normalized 
the light-curve using a time-dependent function with the form 
F baseline = Ao+A\Xt where F baseline is the target flux out of tran- 
sit and t is the time. We also tried higher-degree polynomials and 
logarithmic baselines; this did not improve significantly the fit. 
The linear correction we adopted is not perfect; remaining un- 
certainties in the actual shape of the baseline introduces errors in 
the system parameters derived from the fit (see Sect. 15.1. 41 . The 
baseline and pixel-phase-effect parameters are fitted simultane- 
ously with the transit-related parameters (see Sec. 15.1. 3k this 
allows these effects and their uncertainties to be taken into ac- 
count in the transit parameters determination. The middle panel 
of Fig.|2]shows the light curve after removal of these instrumen- 
tal effects and the lower panel shows the residuals to the fit. 



3. SOPHIE radial velocities 

3.1. Observations 

The first half of the January 2010 transit of HD 80606b 
was visible from Europe so we managed to observe it with 
the SOPHIE spectrograph at the 1.93-m telescope of Haute- 
Provence Observatory in South of France. SOPHIE is a cross- 
dispersed, environmentally stabilized echelle spectrograph dedi- 
cated to high-precision radial velocity measurements (Perruchot 
et al. 120081 Bouchy et al. 120091) . Since the discovery of the 
transit in February 2009, this was the first time that this part 
of the transit was observable from an observatory with high- 
precision spectroscopic capabilities. Unfortunately, the 1.93-m 
telescope started a technical break for maintenance and upgrades 
in November 2009 that was extended up to February 2010. Due 
to the importance of this event, some observations could nev- 



ertheless be performed thanks to the support of the OHP staff, 
wich was mandatory due to the ongoing works on the telescope. 

The observations were secured as part of the second sub- 
program of the SOPHIE Consortium (Bouchy et al. 2009 
Hebrard et al. l2010bl >. The night of the transit, 13 January 2010, 
observations could start just before 23h (UT) after technical is- 
sues were solved and clouds disappeared, and had to be stopped 
4.5 hours later due to cloudy weather. A 5.5-hour sequence of 
good SOPHIE reference observations of HD 80606 could also be 
performed the 15 January 2010 night. 24 and 33 exposures were 
secured during the two nights respectively. The exposure times 
ranged between 5 and 20 minutes, with typical values around 
9 minutes; we tuned it in order to maintain a constant signal- 
to-noise ratio per pixel of ~ 58 at 550 nm despite the weather 
changes (seeing and absorption). 

The measurements were performed with the same setup as 
the one we used for our observation of the February 2009 transit 
(Moutou et al. 2009). We used the fast-read-out-time mode of 
the CCD detector in order to minimize overheads. Observations 
were secured in high-resolution mode allowing the resolving 
power A/AA = 75000 to be reached. The first optical-fiber 
aperture was put on the target and the second one on the sky; 
it allows us to check that no diffuse light was polluting the 
HD 80606 spectra in these Moonless nights. Wavelength cali- 
brations with a thorium lamp were performed with a ~ 2-hour 
frequency each night, allowing the interpolation of the spectral 
drift of SOPHIE for the time of each exposure. A few expo- 
sures were performed with simultaneous thorium-lamp light in 
the second aperture to allow for simultaneous wavelength cali- 
bration; those extra-calibration did not improve significantly the 
radial velocity accuracy. 

3.2. Data reduction 

We used the SOPHIE pipeline (Bouchy et al. |2009l ) to extract the 
spectra from the detector images, to cross-correlate them with a 
G2-type numerical mask, then to fit the cross-correlation func- 
tions (CCFs) by Gaussians to get the radial velocities (Baranne et 
al.[l996 Pepe et al. 2002 ). Each spectrum produces a clear CCF, 
with a 7.82 ± 0.03 km s full width at half maximum and a con- 
trast representing 48.8 + 0.4 % of the continuum. The accuracy 
of the measured radial velocities is typically around 4 m s~ 1 . This 
includes photon noise (typically ~ 2 ms _1 ), wavelength calibra- 
tion (~ 2 m s ), and guiding errors (~ 3 m s -1 ) that produce mo- 
tions of the input image within the fiber (Boisse et al. 2009). We 
also re-reduced the SOPHIE observations presented in Moutou 
et al. (120091 ) in order to have a uniform data set reduced with the 
same version of the pipeline. There was no significant change by 
comparison with the data presented in Moutou et al. (2009), ex- 
cept for one of the reference exposures performed out of the tran- 
sit for which the correction due to the Moon pollution was sig- 
nificantly improved. SOPHIE measurements of HD 80606 per- 
formed in February 2009 and January 2010 have the same prop- 
erties, except the better signal-to-noise ratios for shorter expo- 
sure times obtained in 2010, which is largely due to the primary 
mirror of the telescope that was realuminized in October 2009. 

The SOPHIE radial velocities of HD 80606 are plotted in 
Figs.[3]and|2]together with other data sets. 

In the SOPHIE spectra, the cores of the large Ca n H & K 
absorption lines of HD 80606 at 3934.8 A and 3969.6 A show no 
chromospheric emissions. The level of the Can emission cor- 
responds to logR' HK - -5.3 ±0.1 according to the SOPHIE 
calibration (Boisse et al. in preparation). For a G-type star 
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Fig. 3. Radial velocities and photometry of HD 80606 as func- 
tion of the orbital phase or the time after periastron. Upper 
panel: radial velocities as measured with ELODIE (red open 
squares, Moutou et al. I2009I I, HRS (purple open circles, 
Wittenmyer et al. I20091 . HIRES (Winn et al. 2009a) pre- and 
post-upgrade (green open upward and downward triangles, re- 
spectively) and around the June-2009 transit (green filled down- 
ward triangles), and SOPHIE during the February-2009 and 
January-2010 transits (blue open and filled diamonds, respec- 
tively). Middle panel: same as above, but enlarged around 
the transit phase, lower panel: Spitzer photometry during the 
January-2010 transit. The Spitzer data are binned per 250, i.e. 
one point each 100 sec. The two horizontal dotted-lines show the 
absorption depth expected with the value R p /R* = 0.100 that we 
get and the value R p /R* = 0.103 measured in the visible by Pont 
et al. d20Q9l > and Winn et al. (|2009al >. On the three panels, the 
dashed lines show the models without transit and the solid lines 
show the models with transit (Rossiter-McLaughlin anomaly for 
radial velocities, and absorption feature for the light curve). The 
parameters of the fits are reported in Table Q] The vertical dot- 
ted lines show the mid-transit, the vertical dot-dashed lines show 
the first and fourth contacts, and the vertical dot-dot-dashed lines 
show the second and third contacts. The uncertainty on the tim- 
ing of the mid-transit is 1.5 minute (corresponding to 9 x 10 in 
orbital phase), and about two times larger for the timing of the 
four contacts. 



(B -V - 0.76) with this level of activity, Santos et al. (2000) 
predict a dispersion below 5 ms" 1 for the activity-induced stel- 
lar jitter. According to Noyes et al. (1984) and Mamajek & 
Hillenbrand (2008 ), this level of activity implies a stellar rota- 
tion period P rot > 50 days. 



4. Ground-based photometry 

In order to study its intrinsic variability, HD 80606 was also ob- 
served several nights before and after the January 2010 tran- 
sit. Observations during the night of the transit were prevented 
by bad weather conditions in the observing site, located in 
Portalegre, Portugal. We gathered a total of nearly 23 hours 
of observations in 10 nights spread from December 3, 2009 
to January 25, 2010. The equipment used is a 12-inch f 15.5 
MEADE LX200 telescope, and a SBIG ST8XME 9-micron- 
pixel CCD camera. The field is 28 x 18.7 arcmin 2 and the 
pixel scale is ~ 1.1 arcsec/pixel. Observations were taken 
through a Bessell I filter, using integration times between 40 and 
90 seconds. 

The frames were reduced using standard IRAF routines and 
aperture photometry was obtained for HD 80606, its companion 
HD 80607 and two additional reference stars located in the same 
field, located about 7 and 11.5 arcmin away from HD 80606. 
The size of the photometric aperture was varied on each night, 
depending on seeing conditions, in order to obtain the highest 
possible signal-to-noise ratio. 

In every night, the flux ratio between HD 80606 and 
HD 80607 exhibits a smaller dispersion than the ratio of 
HD 80606 with any of the other two reference stars, or a 
combination of HD 80607 and the additional reference stars. 
Furthermore, the light curves obtained for HD 80607 - using 
the reference stars - also exhibit a larger dispersion than the 
flux ratio between HD 80606 and HD 80607. Therefore, the lim- 
its to the intrinsic variability of HD 80606 were set using only 
HD 80607 as reference star. 

The root mean square for the different nights range from 
about 2.2 to 4.8 mmag, depending mainly on the weather condi- 
tions. During each night, no effect was seen due to variation of 
the sky level, or the position of the stars on the chip or the air- 
mass. On the other hand, the nightly mean of our observations 
exhibits a dispersion of about 1 .7 mmag around its mean value. 

Similar photometric observations of HD 80606 were con- 
ducted the weeks around the January 2010 transit with the CCD 
camera at the 120-cm telescope of Haute-Provence Observatory, 
using the setup used by Moutou et al. (2009). The filter used was 
I'Gmm with a neutral density and we did not defocus. Photometric 
observations started on 6 January and ended on 23 January, with 
a total of eight sequences ranging from 30 to 120 minutes per 
night; in total, 187 images have been acquired and analyzed. 
Exposure times of 10 to 120 seconds have been used to account 
for varying transparency. An observation of the transit during the 
January- 13 night could be performed with this instrument, on a 
coverage similar to this of the SOPHIE observations reported in 
Sect. [3] This ground-based photometry of a part of the transit 
will be studied in a forthcoming paper. 

The flux has been extracted from aperture photometry using 
the GCNTRD and APER procedures (Sect. 12.21) on a 8-pixel radius 
(0.69 arcsec/pixel). The background has been estimated from a 
sky annulus with radii of 10 to 12 pixels. The root mean square 
of flux variations ranges from 1 .0 to 3 . 1 mmag per night, and the 
residual fluctuation has a standard deviation of 3.2 mmag when 
the transit night is excluded. This is significantly larger than what 
was observed in Portalegre, which may indicate a real tendency 
of rotational modulation, or a significant level of systematics. 
No attempt was done to correct for the long-term behavior in 
the data, which appears to be compatible with the expected rota- 
tional period of the star (as discussed in Sect. 16.51 ). 
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We conclude therefore that HD 80606 is photometrically sta- 
ble at the level of a few mmag in the optical range, in the 
timescale of several weeks. 



5. Analysis 

We fitted this whole data set in order to refine the system pa- 
rameters. As a full transit was observed with Spitzer, possible 
systematic effects due to the combination of transit portions se- 
cured with different ground-based instruments are expected to be 
reduced here. In addition, together with the new radial velocities 
secured at phases previously uncovered, the constraints on the 
spin-orbit angle would be better. 

5.1. Combined fit 

5.1.1. Method 

We first performed a combined fit of our Spitzer photometry 
of the January-2010 transit together with the available radial 
velocities of HD 80606. We also included in this combined fit 
the timing constraint on the eclipse as measured from previous 
Spitzer measurements by Laughlin et al. (2009 ) and re-analyzed 
by Gillon (2009a). As those data are not accurate enough to al- 
low the ingress and egress of the eclipse to be significantly mea- 
sured, we used the estimated epochs of the mid-time of these two 
events: HJD = 2 454 424.700 + 0.005 and 2 454 424.775 + 0.005. 

We did not include in this combined fit the radial veloc- 
ity measurements secured during and near the transits. Indeed, 
the Rossiter-McLaughlin observations do constrain the pro- 
jected stellar rotational velocity Vsin/* and the sky-projected 
angle A between the planetary orbital axis and the stellar ro- 
tation axis, but they do not constrain significantly the parame- 
ters that we measure here in the combined fit. The analysis of 
the Rossiter-McLaughlin data is presented below in Sect. 15.21 
and takes into account the results of the combined fit presented 
here in Sect. 15.11 Thus, the combined fit uses the radial veloci- 
ties secured with the instruments ELODIE at OHP (Moutou et 
al. 120091 , HRS at H ET (Wi ttenmyer et al. [2009b , and HIRES 
at Keck (Winn et al. 12009a). This covers a 9.5-year span. The 
SOPHIE data are used only for the Rossiter-McLaughlin fit. We 
note that the ELODIE and SOPHIE radial velocities are absolute 
heliocentric whereas those from HIRES and HRS are relative. 

5.1 .2. Transit light curve of planets on eccentric orbit 

To calculate the transit light curve using a given set of orbital 
parameters (period, orbital inclination, semi major axis in unit 
of stellar radii, eccentricity and longitude of periastron), we cal- 
culated the sky projected distance between the planets and the 
stars center in unit of stellar radius; this last result is then di- 
rectly input to the Mandel & Agol d2002l ) algorithm with limb 
darkening coefficients. From a theoretical model (Kurucz fl979l ). 
the three non-linear limb-darkening coefficients at 4.5 pm as 
defined by Sing (120101 ) has been derived with T e ff = 5500 K 
and logg = 4.5: c 2 = 0.89502981, c 3 = -1.1230710 and 
c\ - 0.46541027. The limb darkening is low at this wavelength; 
Desert et al. (|2009 ) have shown that the uncertainties in the co- 
efficients have no significant effects on the parameters derived 
from the transit light curve. 

It is worth to note that the speed of the light in the system has 
to be taken into account when comparing the times of the tran- 
sit, the eclipse and the periastron. The transit takes place when 



the planet is at ~ 0.29 AU from the star (see Fig. [TJ; this im- 
plies an apparent advance of about 2.5 minutes for the transit. 
Similarly, the planet is at ~ 0.03 AU from the star at the eclipse, 
which implies a delay of 15 seconds. Also because of the planet- 
observer distance decreases during the transit, the speed of light 
correction makes the transit appears to last 8 seconds less than 
its really does. In all the procedure, we implement this speed of 
light correction. 



0.40 f 




-10 -5 5 10 

Time (hours) 



Fig. 4. Plot of the orbital velocity of the planet and orbit direc- 
tion during the transit, a being the angle between the direction of 
the planet motion and the line of sight, sin a is the fraction of the 
planet motion projected onto the sky (upper panel). Although the 
planet velocity significantly decreases during the transit (middle 
panel), the change of the orbital direction has the opposite effect 
such that the planet velocity on the sky is nearly constant (bottom 
panel). Even more, the planet velocity on the sky is nearly sym- 
metrical around the center of the transit, reinforcing the apparent 
symmetry of the transit light curve. On the three panels, the ver- 
tical dotted lines show the mid-transit, the vertical dot-dashed 
lines show the first and fourth contacts, and the vertical dot-dot- 
dashed lines show the second and third contacts. Figure[T]shows 
these parameters drawn on a sketch. 



The orbit of HD 80606b is highly eccentric, and the tran- 
sit takes place ~5.7 days after the periastron when the planet is 
rapidly moving away from the star. During the transit, the star- 
planet distance increases by about 5.8 % between the first and 
the fourth contact, and the orbital velocity decreases by about 
4.1 % (Fig. [4]). With that in mind, a priori, this would suggest 
that the transit light curve could be highly asymmetrical with, 
for instance, an egress lasting longer than the ingress. However, 
during the transit, the direction of the planet motion also varies 
by about the same amount (upper panel of Fig. |4|. If we define a 
as the angle between the direction of the planet motion and the 
line of sight, sin a is the fraction of the planet motion projected 
onto the sky (see insert in Fig.Q]). Between the first contact and 
the last contact, sin a increases by about 4.3 %, exactly compen- 
sating for the decrease of the orbital velocity. As a result, the 
planet velocity on the sky is nearly constant. In addition to the 
low variation of the apparent planet velocity on the sky during 
the transit, it appears that the time variation of this velocity is 
nearly symmetrical around the center of the transit (bottom panel 
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of Fig. |4j». Consequently, the apparent velocity during the egress 
is extremely close to the velocity during the ingress, the later 
being only 0.004 % smaller than the former. This reinforces the 
apparent symmetry of the transit light curve. With the parame- 
ters of the best fit, the ingress duration (time between first and 
second contact) is predicted to last only half a second more than 
the egress duration (time between third and fourth contact) while 
both lasts about 2 hours and 45 minutes. 

This surprising symmetry for the transit light curve of a 
planet on a highly eccentric orbit can be explained by consid- 
ering the projection on the sky of the gravitational force from 
the star to the planet. Because the transited star is the same as 
the attracting star, during the transit the component of the grav- 
itational force projected on the sky is close to zero. As a con- 
sequence, the apparent velocity of the planet projected on the 
sky is nearly constant. Even more, for various positions of the 
planet, the projected force is also symmetrical around the star 
center, explaining the symmetry of the time variation of the pro- 
jected velocity around the center of the transit. As a conclusion, 
in contrary to the common-sense idea, in any configuration the 
transit of a planet in front of its parent star is expected to be 
highly symmetric, even for extremely eccentric orbits (see also 
Winn 2010a). Our Spitzer observations confirm this result in the 
case of the highly eccentric HD 80606b's orbit. 

5.1.3. Free parameters 

Our combined fit of the radial velocities, the Spitzer transit light 
curve and the eclipse timing includes 19 free parameters. We 
list them below, classified in four categories depending on their 
nature and the data set that constrain them: 

1 . two free parameters constrained by photometry only: 

- Rp/R,, the ratio of the planetary and stellar radii; 

- a/R t , the semimajor axis in units of stellar radius; 

2. six free parameters constrained both by radial velocities 
and photometry: 

- a cos i/R*, the "standard" impact parameter, which is dif- 
ferent (for an eccentric orbit) of the actual impact param- 
eter b — d, cos i/R*, where d t is the star-planet distance 
at mid-transit and i is the inclination of the orbit; 

- P, the orbital period of the planet, which is proportional 
to (a//?,) 3 , and to the inverse of the stellar density p» 
from the Kepler third law; 

- 7q, the epoch of the periastron of the planet; 

- e cos to and e sin a>, that constrain the two correlated pa- 
rameters e (the eccentricity of the planetary orbit) and u 
(the longitude of its periastron); 

- K Vl - e 2 , that depends both on the eccentricity e and on 
the semi-amplitude K of the radial velocity variations - 
K actually is constrained by radial velocities only, and 
is proportional to M p and (M,) 2 ' 3 , M p and M, being the 
planetary and stellar masses; 

3. four free parameters constrained by radial velocities only: 

- V„=i_>4, the center-of-mass radial velocities for each of 
the four radial velocity data set used, namely ELODIE 
(Moutou et al. 120091 , HRS (Wittenmyer et al. 120091 , and 
HIRES whose pre- and post-upgrade data are considered 
as two different data sets (Winn et al. 2009a); 

4. seven free parameters linked to the Spitzer light curve ex- 
traction (see Sect. I2.2t : 

- the two parameters Aj = i_>2 for the baseline; 

- the five parameters Kj = \^ for the pixel-phase effect. 



The direct, absolute measurements of M, and R„, and conse- 
quently those of M p and R p , are not feasible from this fit; none of 
those four parameters is measurable independently of the other 
ones from our data set. One exception could be the stellar mass 
that could be directly determined through the semi-major axis a 
thanks the third Kepler law. Indeed, by measuring the delay of 
the eclipse mid-time and the advance of the transit's one that are 
due to the time the light takes to propagate through the HD 80606 
system (~ 2.5 minutes difference between those two events by 
comparison to the ephemeris, see Sect. 15.1.21 ). the semi-major 
axis could theoretically be directly measured in Astronomical 
Units. The uncertainties in the transit and eclipse measured mid- 
times are however of the order of 1 .5 and 4 minutes respectively, 
which implies an accuracy on a of the order of 0.5 AU with this 
approach; this is not constraining here. 

Stellar evolution models remain thus mandatory to esti- 
mate M t . Moutou et al. (120091 1 used isochrones to get M* = 
0.98 + 0.07 M from T eff = 5574 + 50 K, logg = 4.45 ± 0.05 
and [Fe/H] = 0.43 dex. Such accuracy is typical for stellar evolu- 
tion models, from which it is difficult to predict stellar masses at 
better than +10 % (Fernandes & Santos 2004). Using additional 
constraints from the transit light curve, Pont et al. (2009]) and 
Winn et al. (2009a) estimated the stellar mass of HD 80606 to 
0.97 ± 0.04 and 1 .05 ± 0.032 M G , respectively. The 2-cr disagree- 
ment between the two estimates shows that we probably reach 
here the limit of the accuracy that is now achievable. In the fol- 
lowing, we adopt the conservative interval M, = 1.01 ± 0.05 M Q 
that takes into account those two studies. 

5.1 .4. Best parameters and error bars 

We u sed th e Prayer Bead method (Moutou et al. I2004I Gillon 
et al. 120071 ) as applied by Desert et al. (2009) to compute the 
mean values of the free parameters and their statistic and sys- 
tematic uncertainties (see examples in Fig. f5}, together with the 
Levenberg-Marquardt algorithm to provide the best fit at each it- 
eration of our procedure. This method was applied to the Spitzer 
photometry in order to account for possible correlated noise in 
the error budget. We simultaneously applied a bootstrap proce- 
dure to the radial velocity measurements, after having quadrat- 
ically added a systematic uncertainty to the radial velocity data 
sets in order to put to unity their corresponding reduced^ 2 . Thus 
we quadratically added 12 ms _1 to the ELODIE uncertainties, 
and 5.0 and 1.7 ms 4 to the HIRES uncertainties (pre- and post- 
upgrade, respectively). In total, 15000 shifts and fits of transit 
light curves and bootstraps of radial velocity errors were pro- 
duced to derive the set of parameters and to extract their means 
and their corresponding standard deviations. In addition, we per- 
formed additional fits with the light curve considering data start- 
ing at ten different epochs before the ingress; by performing a 
prayer bead on 1500 shifts for each of those ten fits, we could 
estimate the errors caused by the uncertainty in the shape of the 
out-of-transit baseline (see Sect. 12.21) . 

The averaged values we obtained for the seven fitted param- 
eters linked to the Spitzer light curve extraction (Sect. 12.21 ) are: 
Ai = 58818 and A 2 = 6.860 for the baseline, and K\ = 0.2142, 
K 2 = 0.1538, K 3 = 0.1083, K 4 = 0.0729, and K 5 = 0.0050 for 
the pixel-phase effect. However, those parameters could be dif- 
ferent at different epochs and for different pixel locations, so this 
is not clear whether the Aj and Kj values we derived here could 
be applied to other IRAC observations with Waim-Spitzer. 

The best solution of the combined fit is plotted in Fig. [3] with 
a ramp cut up to HJD = -2455 210.325, i.e. the first 1.8 hour of 
the observation was discarded, corresponding to the first 15000 
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Table 1. Parameters for the HD 80606 system 



Parameters 


Values and l-cr error bars 


Unit 


free parameters constrained by photometry only: 




R P /R, 


0.1001 ±0.0006 




a/R, 


97.0 ± 1.6 




free parameters constrained both by photometry and radial velocities: 


a cos i/R, 


1.238 ±0.011 




P 


11 1.4367 ±0.0004 


days 


T Q (periastron) 


2455 204.9 16 ±0.004 


1 1 .11 > 


ecosoi 


0.4774 ±0.0018 




esinw 


-0.8016 ±0.0017 




tfVl -e 2 


171.1 ±0.5 


ms" 1 


free parameters constrained by radial velocities only: 




V ELODIE 


3.7888 ± 0.0023 


km s" 1 


cr(0 - C) 


17.8 


ms" 1 


V HRS 


-0.0193 ± 0.0019 


km s _I 


<r(0 - C) 


6.3 


ms -1 


7 MIRES, pre-ugrade 


-0.1845 ±0.0010 


km s"' 


cr(0 - C) 


5.5 


ms" 1 


'HIKES, post-ugrade 


-0.1827 ±0.0007 


km s _1 


cr{0 - C) 


2.4 




directly derived parameters: 






b 


0.808 ± 0.007 




i 


89.269 ±0.018 


° 


T, (transit mid-time) 


2455 210.6420 ± 0.0010 


HJD 


Transit duration Ti_4 


11.88 ± 0.09 


h 


Transit duration r M = T 3 . 


_ 4 2.78 ±0.10 


h 


e 


0.9330 ± 0.0005 




CO 


300.77 ±0.15 




K 


475.3 ± 2.0 


i 

ms 


P. 


1.39 ±0.07 


g/cm 3 


T e (eclipse mid-time) 


2454424.736 ±0.003$ 


1 1 1 


Eclipse duration ri_ 4 


1.85 ±0.14$ 


n 


derived parameters assuming 


a M* -value: 




M* 


1.01 + 0.05+ 




R* 


1.007 ±0.024 




M p 


4.08 ±0.14 


Mj U p 


R P 


0.981 ±0.023 


Rjup 


P V 


5.4 ± 0.4 


g/cm 3 


a 


0.455 ± 0.008 


AU 


Rossiter-McLaughlin effect: 






Vsin/» 


1.7 ±0.3 


km s _I 


A 


42 ± 8 




V SOPHIE, transit February 09 


3.9162 ±0.0013 


km s _1 


<T(0 - C) 


4.3 


ms" 1 


%/tf£S,post-ugrade, transit June 09 ~ ' 1795 ± °- 001 1 


km s"' 


cr(0 - C) 


0.7 


ms" 1 


* SOPHIE, transit January 10 


3.9018 ±0.0013 


km s"' 


lt(0 - C) 


5.2 


ms" 1 


f: combined value from Pont et al. (20091 and Winn et al. ( 2009a 1 



i : from Laughlin et al. ( 2009 ) 



over the 159232 unbined frames. The upper panel of Fig. [3] 
shows the radial velocities and the lower panel shows the Spitzer 
photometry binned per 250 pixels (the middle panel shows the 
Rossiter-McLaughlin fit that is discussed below in Sect. 15. 2t . 

The derived parameters are reported in Table[T]together with 
their error bars; they are ranked as a function of the way they 
are derived. First are reported the free parameters of the com- 
bined fit that are listed in Sect. 15.1.31 and which are constrained 
by photometry only, then those that are constrained both by ra- 
dial velocities and photometry, and then those constrained by 
radial velocities only. In this last category, the dispersion around 
the obtained radial velocity shifts are also reported. Twelve ad- 
justed parameters of the combined fit are reported here. Then 




0.0985 0.0990 0.0995 0.1000 0.1005 0.1010 0.1015 
Rp/R. 



14- t -(J" + + 




4.910 4.915 4.920 4.925 

Tc_periostre - 2455200 (HJD) 



Fig. 5. Distributions of the parameters obtained for the 15000 fits 
performed with the prayer bead and bootstrap (Sect. 15.1.4b . As 
example, four parameters are plotted here: a/R, vs. R p /R, and P 
vs. JoCperiastron). The mean values and the standard deviations 
are shown. 



Table [TJ shows the parameters that are directly derived from the 
above free parameters of the fit, without any additional hypothe- 
sis. This includes the transit and eclipse timings, the latest being 
obtained from Laughlin et al. (120091 and Gillon (2009a]). The 
following parameters are those that are derived by assuming the 
stellar mass Af» = 1.01 ± 0.05 M (Sect. 15.1.3b together with the 
parameters derived above and through the Kepler third law. The 
last parameter set in Table Q] are those relative to the Rossiter- 
McLaughlin fit that are obtained below in Sect. 15.21 

The dispersion of the Spitzer photometry around the transit 
light curve fit represents 5.3 x 10" 3 of the stellar flux for un- 
binned frames. This is the expected level of the photon noise. 
The amplitude of the correlated noise, as seen for different bin- 
sizes of the light curve, is of the order of 1 .4 x 10" 4 of the stellar 
flux (Fig. |6). A bump in the light curve with an amplitude of 
~ 1 mmag is seen just before the transit mid-time (see Sect l6.1b . 
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It could be instrumental or due to a spot (see Sect. 16. U . We per- 
formed fits without taking into account these points; this did not 
change significantly the derived parameters. 



HD80606 Spitzer/IRAC-ch2 
o.oi oo r 1 1 1 1 1 1 




Fig. 6. Root mean square (RMS) of binned residuals of the 
Spitzer transit light curve as a function of the bin size n. The 
solid, red line is proportional to « as expected if only photon 
noise is considered. 



5.2. Rossiter-McLaughlin fit 

The radial velocity measurements secured during transits have 
been fitted in order to measure the sky -projected angle A between 
the planetary orbital axis and the stellar rotation axis. The data 
we use are the SOPHIE observations of the February-2009 tran- 
sit (Moutou et al. 2009 ) and the new ones we secured in January 
2010 (Sect. E). We also u sed the HIRES data of the June-2009 
transit (Winn et al. 2009a]! . 

We first measured the radial velocity shift of each data set 
by comparison to the center-of-mass radial velocity references 
computed in Sect. 15.11 Radial velocity measurements performed 
near but off the transit are mandatory to constrain those shifts. 
A lack of such reference observations could prohibit an accurate 
measurement of the spin-orbit angle. For example in the case 
of the recent Rossiter-McLaughlin observation of the planetary 
system Kepler-8 (Jenkins et al. 1201 Oi l, the paucity of off-transit 
observations makes difficult to conclude if the apparent asym- 
metry of the Rossiter-McLaughlin shape by comparison to the 
Keplerian curve is due to an actual spin-orbit misalignment or to 
a shift due to another cause, as stellar activity and/or instrumen- 
tal drifts. 

In the case of HD 80606, we used as reference for the 
February-2009 data the nine SOPHIE measurements performed 
the nights before and after the transit night, as well as the nine 
ones performed after the fourth contact the night of the transit. 
For the 2010 SOPHIE data we used the 33 exposures obtained 
the January- 15 night, i.e. two days after the transit. Finally, for 
the HIRES data we used the six June-2009 measurements se- 
cured before and after the night of the transit (the June-2009 
data were excluded from the HIRES post-upgrade data set used 
in Sect. 15. II for the fit of the orbit). The radial velocity shift for 
the three data sets are reported in Table \T\ We found a shift of 
14.4 ± 1.9ms-' between the 2009 and 2010 SOPHIE data sets, 
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Fig. 7. Radial velocities of HD 80606b around transits and their 
Rossiter-McLaughlin fit. The upper panel shows the SOPHIE 
observations around the 14 February 2009 transit (Moutou et 
al. 120091. the middle panel shows the HIRES observations 
around the 5 June 2009 transit (Winn et al. 2009a), and the lower 
panel shows the SOPHIE observations around the 13-14 January 
2010 transit (Sect. [3]). On each of these three panels the data are 
plotted together with the 1-cr error bars. The fits with and with- 
out transit are the solid and dashed lines, respectively. The left 
panels show the data during the nights of the transits, and the 
right panels show all the data secured the nights before and af- 
ter the transits to allow the measurement of the radial velocity 
shift of each data set. The O-C residuals of the fit including the 
Rossiter-McLaughlin anomaly are also plotted. Moreover on all 
panels, the vertical dotted lines show the mid-transit, the verti- 
cal dot-dashed lines show the first and fourth contacts, and the 
vertical dot-dot-dashed lines show the second and third contacts. 



and a shift of 3.2 + 1.3 ms" 1 between the HIRES measured in 
June 2009 and the other post-upgrade HIRES data secured more 
than one year earlier. These shifts could result from a combina- 
tion of instrumental effects, unknown component in the system 
and/or activity-induced jitter. 

The reference observations secured near each of the three 
transits are plotted in the right panels of Fig. [7] Using these off- 
transit radial velocities, we also computed the uncertainties we 
had to add to the radial velocities tabulated error bars in order 
to put to unity the reduced x 2 corresponding to the Keplerian 
fit. We thus quadratically added 1.0 ms" 1 and 2.5 ms -1 to the 
SOPHIE measurements secured in 2009 and 2010 respectively, 
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and 1.7 ms -1 to the HIRES measurements (in agreement with 
the other HIRES post-upgrade data used above, see Sect. 15.1.4b . 

We model the Rossiter-McLaughlin anomaly shape using the 
analytical approach developed by Ohta et al. ((2005 >. The com- 
plete model has 14 parameters: the stellar limb-darkening linear 
coefficient e, the transit parameters R p /R», a/R, and i, the orbital 
parameters (P, To, e, u, K), the three radial velocity shifts mea- 
sured above, and finally V sin/* and A. We computed e = 0.722 
in the wavelength range 5300-6300 A using the same method as 
in Sect. 15.1.21 The transit and orbital parameters are fixed from 
the results obtained in Sect. 15.11 their uncertainties are negli- 
gible by comparison to those of the two main free parameters 
of the Rossiter-McLaughlin fit: A, which is constrained by the 
asymmetry of the anomaly, and V sin /*, which is constrained by 
its amplitude. 
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Fig. 8. Confidence interval contours for the A and V sin /* values 
from the Rossiter-McLaughlin fit. 



The best fit is plotted in the middle panel of Fig. [3] and in 
Fig. It is obtained for A = 42° and Vsin/* = 1.7kms _1 ; its 
X 2 is 1 17.8 for 1 17 degrees of freedom. The confidence interval 
contours estimated from^ 2 variations (Hebrard et al. 2002J for 
the two correlated A and V sin I* parameters are plotted in Fig. [8] 
The uncertainties obtained this way are ±6° and ±0.2 km s _1 , re- 
spectively. We increased them in order to take into account for 
the uncertainties in the radial velocity shifts of the three data sets 
with respect to the Keplerian curve computed in Sect. 15.11 Our 
final values are A = 42° ± 8° and V sin/, = 1.7 ± 0.3 km s" 1 . 
We checked that the uncertainties on the parameters derived in 
Sect. 15.11 and used here in the Rossiter-McLaughlin fit imply 
negligible uncertainties on A and Vsin/*; similarly, the uncer- 
tainty in the e limb-darkening coefficient is negligible. 

The adopted values are reported in Table [1] together with 
the dispersion of the radial velocities performed during the tran- 
sit with respect to the Rossiter-McLaughlin fit. The dispersion 
is in the range 4 - 5ms~' for SOPHIE data; we have a simi- 
lar dispersion around the Keplerian curve for the SOPHIE ref- 
erence measurements performed before and after the transit. For 
the HIRES data, the dispersion around the Rossiter-McLaughlin 
fit is 0.8 ms -1 , whereas it is two times larger (1.7ms _1 ) for the 
June-2009 HIRES measurements performed before and after the 
transit. So it could remain a slight systematic in the HIRES ra- 
dial velocities secured during the transit. However this effect is 



small, and is detected here due to the particularly high accuracy 
of these measurements. 



6. Discussion 

6.1. Warm-Spitzer transit light curve 

This planetary transit observations is among the first ones se- 
cured with Spitzer in its post-cryogenic mission. This shows that 
accurate transit light curves can be obtained in this second part 
of the observatory mission despite enhanced ramp and pixel- 
phase effects. Our Waim-Spitzer transit light curve of HD 80606 
has lower uncertainties than those obtained from the ground. The 
typical durations of the pixel-phase effect (70 minutes) and the 
planetary transit (12 hours) are on time-scales different enough 
to avoid significant uncertainties on the derived system parame- 
ters due to this instrumental effect. The transit light curve is well 
detected, without extra signatures of transiting rings or satellites. 
According to our accuracy, the signature of an hypothetic satel- 
lite would have been detected up to a magnitude depth slightly 
below 1 mmag; this corresponds to an upper limit of 2 R e on 
the radius. 

A bump in the light curve is seen just before the transit mid- 
time with an amplitude of ~ 1 mmag. It could corresponds to the 
planet occulting a dark spot during the transit (Pont et al. 2007 
Rabus et al. 120091 ). HD 80606 is not an active star, so the prob- 
ability of such event is low. The stability at the level of a few 
mmag of the stellar flux observed in the optical (Sect. HI would 
correspond to a stability below 1 mmag at 4.5 /mi, assuming 
dark spots are 1000-K cooler than the stellar surface (Desert et 
al. 12010b ; so the optical photometry cannot exclude that the fea- 
ture detected in the Spitzer light curve is due to a spot. One can 
note that at the time of this feature (near 137 hours after the 
periastron in Fig. [3), a simultaneous feature is also seen in the 
SOPHIE radial velocity data. This could argue in favor of the 
interpretation of this event in term of phenomenon at the sur- 
face of HD 80606. This seems however unlikely, since a photo- 
metric feature with this low flux amplitude should a priori not 
produce such high radial velocity effect. The radial velocity fea- 
ture is more likely an instrumental systematic. The feature in the 
Spitzer light curve is detected at a maximal value of the pixel- 
phase effect (Fig. |2j so we should remain cautious about instru- 
mental effects. We monitored the coordinates of the target on the 
subarray through the observation in order to test if the bump cor- 
responds to a particular area of the detector. This is not the case: 
at the epoch of the bump, the target is located on a position of 
the detector where the target passes numerous times before and 
after this event, and where the pixel-phase effect apparently is 
well corrected by our procedure. So we found no strong reasons 
to particularly favor an instrumental effect to explain the pres- 
ence of this feature in the Spitzer light curve. The most likely 
explanation is the presence of a small star spot on the surface 
HD 80606, above which the planet is transiting. 

Our data set results from the first joined campaign during 
which space-based photometry and high-precision radial veloc- 
ities were carried out simultaneously. If large stellar spots are 
detected during a similar campaign, this would help to better 
understand the effect of stellar activity on radial velocity mea- 
surements. Stellar jitter indeed introduces noise in measured ra- 
dial velocities and is a significant limitation in high-precision 
velocimetry (see, e.g., Saar & Donahue 1997, Santos et al. 2000 
Boisse et al. 120091 Queloz et al. 120091 Hebrard et al. l2009ab . No 
stellar spots large enough to be spectroscopically detected were 
apparently seen here. 
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6.2. Comparison with previous measurements 

The HD 80606's system parameters that we report in Table Q] 
have better accuracy by comparison with previous studies 
( Moutou et al. [20091 Pont e t al. [20091 GiUon I2009al Winn et 
al. I2009al Hidas et al. I20101 >. With respect to the groun d obser - 
vation of nearly all the transit phases used by Winn et al. ( 2009a ), 
the uncertainties presented here are better by factors two to five. 
Exceptions are the parameters e, cj and K for which the error 
bars are not significantly reduced. This is due to the fact that 
most of the constraints on these three parameters come from the 
radial velocity on a 9.5-year time span and already used by pre- 
vious studies. 

For most of the parameters, the revised, more accurate val- 
ues are in agreement within 1 cr with those previously published. 
There is a 2 - cr disagreement on the inclination i of the orbit, 
which is found in our study to be slightly lower; this implies a 
more grazing transit, a larger impact parameter b, and a slightly 
longer duration for ingress and egress (~ 10 minutes longer). 
The a//?*-ratio is also found to be 5 % smaller, a shift by 2cr 
according the error bars from Winn et al. d2009al l. 

The projected stellar rotational velocity we found, V sin I* = 
1.7 ± 0.3 km s" 1 , is slightly larger than the value obtained by 
Winn et al. (2009a) from the Rossiter-McLaughlin fit of their 
HIRES data: Vsin/» = 1.12+9-^ km s _1 . Whereas we also used 
their HIRES in our analysis, our different result is explained 
by the fact that we use additional SOPHIE data and we mea- 
sured a smaller R p /R* ratio (see Sect. 16.3b . In addition, a dif- 
ference with the study by Winn et al. (2009a) is the shift of 
3.2+ 1.3 ms~' that we found between the HIRES data secured 
near the transit and the other ones (Sect. 15.2b ; this implies a 
slightly larger amplitude for the Rossiter-McLaughlin anomaly, 
thus a larger V sin/». The value V sin/* = 1.7 ± 0.3kms~' that 
we measured nonetheless agrees with those found from syn- 
thetic spectral fitting: V sin I„ = 1.8 ± 0.5 and 2.0 + 0.5 km s' 1 
by Valenti & Fischer (120051 ) and Winn et al. (2009a), respec- 
tively. The discrepancy noted and discussed by several authors 
(Winn etal|2005l Triaud et al. l2009l Hirano et al. l20T0l Simpson 
et al. 120101 ) between the V sin /» measured from the Rossiter- 
McLaughlin effect and from the spectral modeling of line broad- 
ening is negligible for such slow-rotating star like HD 80606. 

The stellar radius we measured, — 1 .007 + 0.024 R , is 
slightly larger than this obtained by Winn et al. (2009a), namely 
R+ = 0.968±0.028 R Q . Both values remain in agreement with the 
one deduced by Moutou et al. (2009 ) from relationships between 
stellar radius, luminosity, temperature, gravity and mass, namely 
R» =0.98 + 0.07 R . 

The radius ratio R p /R* and the timing of the transit that we 
found significantly differ from the values derived from ground- 
based observation; this is discussed below. 

6.3. Planetary radius 

The radius ratio we found from the transit observed with Spitzer 
is R p /R t = 0.1001 + 0.0006. By comparison, Winn et al. (I2009at 
measured R p /R, = 0.1033 + 0.0011; Pont et al. ( [20091 found 
the same radius ratio than Winn et al. (2009a) but with an un- 
certainty three time larger, mainly due to the lack of ingress ob- 
servation. The radius ratio measured in the infrared from space 
is thus ~ 3% smaller than the one measured in the optical from 
ground. This is a 3 - cr difference according the error bar re- 
ported by Winn et al. d2009al ). The two horizontal dotted lines 
in the lower panel of Fig.[3]show the absorption depths expected 



with the two values; the infrared radius ratio clearly is smaller 
than the optical one. 

The uncertainties in the background measurement (Sect. 1221 
are too small to account for such a radius difference. A varia- 
tion of the stellar brightness due to spots could explain time- 
variations in the measured radius ratio (Desert et al. 2010). The 
upper limit of the brightness variations of the non-active star 
HD 80606 (Sect. [4| is however too small to explain a ~ 3 % 
radius-ratio variation. The planetary thermal emission could pro- 
duce an underestimation of the measured radius ratio (Kipping 
& Tinetti 2009). The emission level measured by Laughlin et 
al. ([2009) at 8 yum near the periastron is larger than the planetary 
emission expected at 4.5 fim at the transit; even with this overes- 
timated planetary flux, the effect would be negligible according 
the error bars on the radius ratio we measure at 4.5 fim. 

The question of possible interpretation in terms of differen- 
tial atmospheric absorption could be raised. The planet appar- 
ent radius as a function of wavelength has been shown to fol- 
low the equation dR,,/dA = Hd\ncr/dA (Lecavelier des Etangs 
et al. 2008a, 2008b), where A is the wavelength, cr is the cross 
section of the mean atmospheric absorber, H = kT/fig is the 
atmospheric scale height, k is the Boltzmann constant, T is the 
temperature, fi is the molecular mass, and g is the planet grav- 
ity which is about 100 ms~ 2 . If the variation of radius is due to 
variation of absorption by haze in the atmosphere as in the case 
of the HD 189733b (Pont et al. 120081 , even assuming Rayleigh 
scattering which produces the steepest variation of absorption as 
a function of wavelength, a temperature of about 5 000 K is re- 
quired to interpret the present measurements, much higher than 
actually measured (Laughlin et al. 2009 ). Indeed, because of the 
large planetary mass, at a typical temperature of 1000 K the scale 
height is only 60 km. The difference of measured radius is about 
40 times this scale height, and is therefore unlikely due to atmo- 
spheric differential absorption. 

We found no astrophysical interpretations able to explain the 
radius difference between optical and infrared wavelengths. It 
is more likely that the error bars are slightly underestimated, 
maybe in the ground-based composite light curve. We note that 
as we obtain a larger stellar radius than Winn et al. (20093), we 
obtain a similar planetary radius despite the different radius ra- 
tio: R p = 0.981 + 0.023 R Jup . 

6.4. Transit timing 

The mid- times of the transit and the eclipse given in Table [T] are 
those that are measured. This means that the transit mid-time 
T t we report is in significant advance by about 2.5 minutes in 
comparison to this predicted from the epoch Tq we report for 
the periastron of the planet (see Sect. 15.1.21 ). Indeed, this latest 
time is in the referential frame of the radial velocities, which are 
these of the star. Similarly, the epoch T e of the eclipse reported 
in Table[T]is delayed by about 15 seconds by comparison to the 
ephemeris computed from To 

The mid-time of the transit we obtain is accurate at the level 
of 1 .5 minute, whereas the accuracy is 6 minutes for the epoch of 
the periastron. Laughlin et al. (2009) obtained a 4-minute accu- 
racy for the mid-time of the eclipse. These accuracies are high, 
especially when they are compared with the long duration of the 
transit and the long orbital period. Using radial velocities alone, 
even if they span a time as long as 9.5 years, the error bars on 
the eclipse and periastron timing are three times larger. For the 
transit timing, the radial velocities can not predict it at better 
than two to three hours. Transit and eclipse detections allow here 
more accurate timings. 
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This high-accuracy offers opportunity to look for possible 
transit timing variations (TTVs). According to the mid-transit 
time and the period we found, we obtain for the February and 
June 2009 transits a mid-time that is ~ 23 minutes earlier than 
the T, time measured by Winn et al. (2009a) from these events. 
Their accuracy on this timing was ±7 minutes, so the disagree- 
ment is at the level of 3 cr. If it is not caused by underestimated 
systematic uncertainties, such a difference in the transit timing 
could in principle be due to the presence of additional bodies 
in the system, such as satellites of the transiting planet or addi- 
tional planetary-mass bodies in the system (see, e.g., Holman & 
Murray [20031 Agol et al. l2005l Nesvorny & Beauge l20T0l i. 

We explored a small region of the parameter space of a hypo- 
thetical additional planet in search for a few examples that could 
explain a ~ 20-minute difference between two nearby transits. 
For this, we performed a series of 3-body simulations by inte- 
grating the equations of motion using the Burlisch-Stoer algo- 
rithm implemented in the Mercury6 package (Chambers 1999 ). 

Planets in resonant orbits could explain such TTV even with 
masses low enough to prohibit their detection with the available 
radial velocities. For example, a 15-Earth-mass planet in a circu- 
lar 4:1 -resonant orbit produces the adequate timing anomalies, 
with many pairs of transits exhibiting a ~ 20-minute discrep- 
ancy. The simulated timing variations exhibit an amplitude of 
80 minutes in this case. A 0.17-Jupiter-mass planet at the 6:1 
resonance produces similar results, although the overall ampli- 
tude is reduced to ~ 40 minutes. Both cases would imply radial 
velocity variations with semi-amplitudes K ^ 1.5 and 4ms~ 1 , 
respectively, so we cannot exclude those two cases: the first one 
would be undetectable in the available radial velocity data set, 
whereas the second one would be at the limit of detection. More 
massive planets on non-resonant orbits could be excluded. For 
example, a planet located in a circular orbit at 3 AU and hav- 
ing a mass of 1 .7 Jupiter masses would imply a small number of 
pairs of nearby transits exhibit a variation of the order of 20 min- 
utes. It would imply radial velocity variations of K ^ 30 ms -1 
which is detectable with the available radial velocities. In all 
these cases, the short-term variations (in timescales of the order 
of 2000 days) exhibited by the orbital parameters of the transit- 
ing planet remain below the level of precision reached from our 
combined fit (Sect. I5.lt . However, the variations in timescales 
of the order of 10 5 ~ 6 days are substantially larger in all cases. In 
particular, some 4: 1 -resonant cases could be unstable after a few 
thousands of years. 

We also explored the possibility that the observed discrep- 
ancy were produced by a satellite to the transiting planet. In the 
more favorable case of having a satellite orbiting the planet at 
the Hill's radius, we found that its mass must be around 3/100 
that of HD 80606b for it to be the cause of a ~ 20-minute delay 
in the transit occurrence, corresponding to 40 Earth masses. 

6.5. Spin-orbit misalignment 

We confirm the spin-orbit misalignment in the HD 80606 system 
and reduce the uncertainty in its measurement: A = 42 + 8°. 
This allows scenarios with aligned spin-orbit to be rejected with 
a high-level of confidence. Such misalignment would imply an 
asymmetry in the Rossiter-McLaughlin anomaly as the star is 
not a perfect sphere but is slightly elongated at the equator due to 
its rotation. However, the accuracy of the data is not high enough 
to allow such tiny effect to be detected. 

The A angle measures the sky-projected angle between the 
planetary orbital axis and the stellar rotation axis. Its actual value 
remains unknown, as the inclination 7, of the stellar rotation axis 



is undetermined. In cases where A = is measured, it is rea- 
sonable to assume /* ^ 90° and an actual spin-orbit alignment. 
In cases as HD 80606 where A is significantly different from 0, 
there is certainly no reasons to assume the stellar rotation axis is 
parallel to the sky plan. 

The projected stellar rotational velocity we get from the 
Rossiter-McLaughlin fit is V sin 7* = 1. 7 + 0.3 km s -1 . According 
to our measured radius for HD 80606, /?* = 1.007 ± 0.024 R , 
and to the relation from Bouchy et al. (2005), this V sin/» trans- 
lates into a stellar rotation period P Iot / sin/, = 30 + 5 days. So 
the stellar rotation period apparently is shorter than 40 days. This 
is shorter than the limit P I0t > 50 days we get in Sect. l3.2l from 
the low activity of HD 80606. This estimation based on the ac- 
tivity, however, is not accurate especially for long rotation pe- 
riods. These two values suggest a rotation period in the range 
40-50 days, so a /*-angle near 90°. The value A = 42+8° that we 
measure is thus probably close to the value of the misalignment 
angle without projection effect. So the actual angle between the 
spin-axis of HD 80606 and the normal to the planetary orbital 
plane is about 40°. 




J WASP- 17 



5 10 15 20 25 

Planetary mass (Jupiter masses) 

Fig. 9. Sky-projected /l-angle between the planetary orbital axis 
and the stellar rotation axis as a function of the planetary mass, 
for 21 published systems (see references in Sect. 16. 5t . Systems 
with |A| > 30° are in filled, red symbols. HD 80606 is marked by 
a circle. 



The first case of a planetary system with a stellar spin mis- 
aligned with the normal of the planetary orbit was discovered 
by Hebrard et al. (2008) in the XO-3 system. This result was 
hereafter confirmed by Winn et al. (2009b), who however found 
a lower /l-value. HD 80606 was the second system reported to 
have a spin-orbit misalignement by Moutou et al. ([2009). This 
result was subsequently confirmed and refined successively by 
Pont et al. d2009l . Gillon d2009at . Winn et al. d2009at and finally 
by the present study. Other planetary systems with significantly 
misaligned spin-orbit have been reported since then: WASP- 14 
(Johnson et al 120091 . WASP- 17 ( Anders on et al. l20T0b . HAT-P- 
7 ( Winn et al. l2009cl Narita et al. l2009at . and CoRoT-1 (Pont et 
al. [2010). There are thus now six known misaligned systems. 
Three other systems may be misaligned, but the large uncer- 
tainties in the reported /l-values prohibit definitive conclusion: 
WASP-3 (Simpson et al. l20T0l), TrES-1 (Narita et al. [2007T) . 
and CoRoT-3 (Triaud et al. 2009). The recent case of Kepler- 
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8 (Jenkins et al. 1201 Oi l is presented as moderately misaligned 
but this requires confirmation (see Sect. 15.21) . On the other 
hand, eleven systems are apparently aligned, namely HD 209458 
(Quelozet al. HEM Witte nmyer et al. l2005l HD 189733 (Winn 
et al. 120061 Triaud et al. |20U9l Collier Cameron et al. 2009), 
HD 149026 (Wolf et al. [20071 Winn & Joh nson in prep.), 
H D 17156 (Cochran et al. [20 081 Ba rbieri et al. [20091 Na rita et 
al. |2009bl ), H AT-P-1 (Johnson l2008b . HAT-P-2 (Winne t al. l2007l 
Loeillet et a l. 120081 HAT-P-1 3 (Winn e t al. 1201 Obi ). WASP-6 
(Gill on et a l. [2009b ), TrES-2 (Winn et al . 120081 TrES-4 (Narita 
et al. l20"T0l ). and CoRo T-2 (B ouchy et al. l2008b . 

Fabry cky & Winn (2009 ) have shown that the whole sample 
of the measured /i-angles could be well reproduced with a bi- 
modal distribution, by assuming that a fraction of the orbits have 
random orientations relative to the stars and the remaining ones 
are perfectly aligned. The systems of misaligned orbits would 
be those which experienced gravitational interactions between 
planets and/or stars, such as Kozai migration. 

In Fig.[9]we show the measured /l-angles as a function of the 
planetary masses for the 21 systems with published measure- 
ments reported above. This plot suggests a scenario with three 
distinct populations. Indeed, for the planets with masses similar 
to that of Jupiter, most of the spins are aligned with the orbits; 
this is expected for planets that formed in a protoplanetary disk 
far from the star and that slowly migrated closer-in at a later 
time. A small fraction of these Jupiter-mass planets however ex- 
hibits large /l-angles; this could be the signature of the second 
population in the scenario by Fabrycky & Winn (2009), these 
which experienced gravitational interactions and that apparently 
are less frequent. These planets also seem to be those than can 
have extreme /1-values. 

A third population could be formed by the large-mass plan- 
ets. Indeed, most of them exhibit a misaligned spin-orbit (see 
also Johnson et al. 2009), maybe all of therrQ. This is a priori sur- 
prising as one can expect the most massive planets are the ones 
for which exciting the inclination is more difficult. If indeed such 
misalignments are frequent for the high-mass-planet population, 
this suggests a different evolution scenario for them. Maybe the 
more massive planets could not really slowly migrate because 
of the interactions with the disk. In that case, only more severe 
interactions with another planet or a star could be the cause of 
migration for massive planets, such interactions also affecting 
the inclination of the orbit. We also note that the /l-angles seem 
lower for the massive planets than those of the misaligned plan- 
ets with lower masses; this suggests as well a different scenario 
for the processes able to modify the inclination of low- and large- 
mass planets. 

Rossiter-McLaughlin observations of other systems with 
transiting massive planets should be performed to confirm or not 
that there are preferentially tilted. 

7. Conclusion 

We presented an observation of the 12-hour-long transit of the 
highly eccentric, 111 .4-day -period exoplanet HD 80606b per- 
formed in January 2010. The transit light curve we present is 
among the first ones carried with the post-cryogenic Spitzer. Its 
shows systematic effects stronger than those seen in the Cold- 
Spitzer but the accuracy remains clearly better that ground- 
based observations. Together with the SOPHIE measurements 
acquired at the same epoch, this is one of the first observational 

3 The only exception, HATP-P-2, could actually be also misaligned 
as the inclination of the stellar rotation axis is unknown. 



campaigns performed simultaneously in radial velocities and 
space-based high-accuracy photometry. With previously avail- 
able data sets, this allows the parameters of this system to be 
significantly refined thanks to combined fits. There is a possi- 
ble detection of a variation in the transit timing, which has to 
be confirmed by additional observations of HD 80606b transits, 
from ground or space. A dark spot was also possibly detected 
on the surface of this inactive star. The spin-orbit misalignment 
is clearly confirmed and the A angle is accurately measured. As 
most of the massive planets for which this angle is measured, the 
orbit of HD 80606b is misaligned with the equatorial plan of its 
host star. This suggests a separate evolution scenario for massive 
planets in comparison with Jupiter-mass planets. 
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